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The nature of dynamics of opinion formation modeled as a decision-by- 
majority process in complex networks is investigated using eigenmode analy- 
sis. Hamiltonian of the system is defined, and estimated by eigenvectors of the 
adjacency matrix constructed from several network models. The eigenmodes 
of initial and final state of the dynamics are analyzed by numerical studies. 
We show that the magnitude of the largest eigenvector at the initial states 
are key determinant for the resulting dynamics. 

1 Introduction 

Many phenomena related to spreading, synchronization and collective dynam- 
ics have been studied from the viewpoint of complex networks. These studies 
have revealed some significant influences of network structure on the phenom- 
ena, and contributed to the development of complex network analysis [1]. In 
such studies, a simple interaction model has been utilized. This is because 
complex phenomena often emerges from simple interactions. 

We focus on dynamics of opinion formation in networks. One of the ap- 
proaches for this dynamics is referred to spin systems. In such systems, local 
interactions are modeled using spin-like variables located at the vertices of 
the networks. Complex dynamics have been explored both analytically and 
computationally. Strong dependency on properties of the networks has been 
pointed out [2, 4, 6, 7]. In some networks, the final consensus state does not 
always have all the vertices with the same opinion [4, 6]. Moreover, it has 
been found that history of convergence to a steady state or a quasi-steady 
state depends on the topology of networks, as well as the local rule [5]. We 
also studied an effect of initial conditions on Glauber dynamics for Ising model 
[9] . We found that non-random initial conditions play a key role to determine 
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the final consensus state, and that the final fraction of each opinion r' is a 
function of initial fraction of the corresponding opinion r. 

In this paper, we will manifest the mechanism that all the vertices with 
the same opinion does not always appear in the final state, and consider the 
origin of the specific r — r' relations under different initial conditions. We 
assume that the mechanism is related to the eigenmodes correspond to each 
eigenvectors of adjacency matrix of networks. Dynamics of opinion formation 
is analyzed using eigenmode analysis. 



2 Models 

Let Vi denote a vertex in a network. For each vertex, there are two possible 
states which represent two opposite opinions. These states are represented 
by a spin- like variable: <Ji(t) = ±1 at a time t. Oi{t + 1) is determined only 
by the values of neighboring vertices of Uj. We adopt a local rule driven by 
decision-by-majority process. It is described by Eqn. (1) 

<7i(t + 1) = sgn|^Oija-j(i)| , (1) 

where a>ij is a component of the adjacency matrix A of the network, which 
takes the value of 1 if an edge exists between vertices Vi and Vj, otherwise 0. 
<7i(t) is updated synchronously at each step t; the values at all vertices are 
updated simultaneously as t progresses. In our previous studies [8, 9], we found 
that depending on initial conditions and network structures, the state of the 
system at t — oo from Eqn. (1) becomes either a fully ordered state in which 
all vertices have the same state, or a metastable state in which two states 
co-exist. In this paper, we analyze the final states by eigenmode analysis. 

We assume that the system controlled by the local rule (Eqn. (1)) evolves 
as the energy becomes lower, and the most stable state appears. First, we 
define a Hamiltonian H of the system as follows: 

H = --^aiajOij, (2) 

i,3 

where Gi and aj are the states of vertices Vi and vj, J is a positive constant. 
This Hamiltonian is regarded as the energy of our system. Equation (2) means 
that the local energy becomes lower if vertices with the same state are adja- 
cent. Secondly, we introduce another expression of H using vector form. Let 
s be the spin state vector at time t; s — {ci(t) 02 (i) • ■ ■ cr n {t)) , where n 
is the number of vertices. The spin state vector s is expanded in the vector 
subspace spanned by eigenvectors of A as s = Y^i=i ° iVi ^ wnere v % is the ith 
eigenvector of A, and Ci is the ith coefficient corresponding to the eigenvector. 
After some calculus, we obtain 
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(3) 



where Aj is the ith eigenvalue of A. For simplicity, we deal with undirected 
and unweighted networks, in which A becomes a real symmetric matrix, and 
all the eigenvalues are real. Thus, the eigenvalues are labeled in descending 
order without loss of generality; Ai > A 2 > • • • > A„. Considering that Ai 
is the largest eigenvalue, it can be expected that the most stable state is 
excited by such the s that gives the maximal C\. It can be shown that the 
eigenmode corresponds to the largest eigenvalue dominantly appears in the 
system. However, we cannot explain a metastable state (in which two states 
co-exist) by this eigenmode. In the following section, we numerically study 
the eigenmodes using several models of complex networks. 

3 Numerical Studies 

3.1 Network Models and Initial Conditions 

In this paper the following network models are studied: Erdos-Renyi ran- 
dom graph (ER), Watts-Strogatz (WS) model Barabasi- Albert (BA) model, 
Klemm-Egmlz (KE) model and Connecting Nearest Neighbor (CNN) model. 
The number of vertices is n = 3000, and (k) = 10 for the average degree. Each 
model has different structural characteristics. See Ref. [1] for the details. In 
such complex networks, vertices are not interconnected homogeneously. We as- 
sume that the resultant dynamics depends on the distribution of initial state 
denoted by Sq owing to this heterogeneity. In order to study the dependency 
of so, we consider arbitrary distributions of initial states according to several 
types of centrality; degree centrality and closeness centrality of the vertices, 
as well as random distribution. See Ref. [3] for the detail of these centrality 
measures. At t = 0, fraction r and distribution of the vertices which a takes 
+ 1 are determined. Corresponding to the centrality measures, the rn vertices 
with the largest centrality are assigned as cr(0) = +1, while the remaining 
(1 — r)n vertices are assigned cr(0) = — 1. 

3.2 Relationships between Eigenmodes and Dynamics 

First, we consider the initial state given by a random distribution. The square 
of the first coefficient in the eigenmode expansion of s normalized by n {c\/n), 
which is derived from the largest eigenvector and minimizes the Hamiltonian 
is examined. The solid lines plotted in Fig. 1(A) and (B) present c\/n versus 
the initial fraction of +1 for Sq and Sqo, respectively, where denotes the 
final state. Fig. 1(C) shows the fraction of +1 in the final state. In the initial 
state, the largest eigenmode is prominent, except in the region of r ~ 0.5. In 
that region, c\/n becomes close to zero. We have checked that the eigenmodes 



4 Makoto Uchida and Susumu Shirayama 



associated with smaller eigenvectors are dominant. Note that J2i c i/ n = 1- 
We also find the differences of the magnitude of c\/n among the networks. In 
the ER, WS and BA networks, it is close to one, while it is about 0.1 for the 
KE and CNN networks. At t = oo, c\/n is superior to the other associated 
with smaller eigenvectors in most range of r, as shown in Fig. 1(B). The 
distributions of magnitude of c\jn among the networks are different. In the 
ER and the WS networks, c\jn is close to one. We have checked that the 
magnitude of c^/n for the KE network is almost same order of the c\/n. From 
Fig. 1(B) and (C), relationships between r — c^/n and r — r' relations can be 
exploited according to the change caused by the variation of r. In the range 
of r where the largest eigenvector is dominant, all vertices are in the stable 
state, that is, they have a single opinion. The range that the two opinions 
coexist in the final state can be found from r — r' relations. In such range of 
r, c\/n becomes smaller than the magnitude at the other range. 

Secondly, we examine the initial state distribution determined by degree 
centrality and closeness centrality, whose results are also represented in Fig. 1. 
The range of r where c\/n becomes smaller varies according to the initial 
conditions and network structures. Let r c be the value of r which c\/n takes 
nearly zero. In the case of the ER and WS networks, r c is slightly less than 
0.5. The value of r c for the other networks is much less than 0.5. In Soo, cf/n 
is dominant in most range of r, while the value of c\/n is relatively small 
compared to the ER and WS networks. The range of r where c\/n ~ in 
agree with that in Sq. Moreover, r' in clearly describe the characteristics: 
Around the value of r c in Sqc, r' rapidly increases as r. Furthermore, in the 
range of r where the dynamics ends up in a metastable state with coexisting 
two opinions, c\/n has much smaller value, and it is assumed that the others 
associated with smaller eigenvectors have relatively large value. 

4 Discussion and Conclusion 

The eigenmodc of the largest eigenvector implies the mode that all vertices 
have the same opinion. This corresponds to the initial state that the fraction 
of two opinions r is r = or r = 1 in s . On the other hand, one can 
consider a particular r where s and the largest eigenvector v\ are almost 
orthogonal. Since the components of i>i are non-negative and the components 
of s is ±1, s and v\ becomes almost orthogonal if randomly selected half 
of the components of Sq are +1, while the others are — 1. This is because the 
eigenmode of the largest eigenvector cf/n becomes c\/n ~ in s at r ~ 0.5, 
if the initial state is randomly distributed. 

Contrary to this, if the components of Sq whose corresponding components 
of v\ are large, are preferentially assigned +1 at the initial state, the inner 
product ci of s and v\ becomes c\ ~ although +1 is less than half in 
the components of So, that is, r < 0.5. The value of r is less than 0.5 when 
the initial states are distributed according to the centralities. This suggests 
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that a certain correlation exists between the values of the elements of i>i and 
such centralities. In the final state Soo the fractions of two opinions r' steeply 
increase around the values of r where c\/n ~ in s . Such r agrees with the 
value where s and v\ are orthogonal. This fact result suggests that the final 
states and the values of r at which r' have rapid transition can be estimated 
by analyzing the magnitudes of the largest eigenmodc on sq. 

In summary, we have performed an eigenmode analysis on dynamics of 
the decision-by-majority process on complex networks. Relationships between 
initial conditions and final states of the dynamics have been analyzed by 
eigenvectors of adjacency matrices. First, we confirmed analytically that the 
stable state in which all vertices have the same opinion corresponds to the 
first eigenmodc of the system. Then, we analyzed the relationships between 
the dynamics induced by arbitrary initial conditions and eigenmodes of net- 
works computationally. It has been shown that, the values of r at which the 
magnitude of the largest eigenvector becomes extinct in the initial state, gives 
r' a rapid transition in the final state. If the magnitude of the second or 
less eigenvector is stronger than the largest eigenvector in the initial state, a 
metastable state arise with two opinions coexisting at the final state. From 
this fact, the final state of the dynamics can be estimated by the magnitudes 
of the largest eigenmode at the initial state. In conclusion, we have thus shown 
that the eigenmodc analysis gives us a clue for understanding and predicting 
an evolution of dynamics on complex networks. 
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Fig. 1. Numerical results of (A) Normalized magnitudes of the eigenmode of the 
largest eigenvectors on so, (B) Normalized magnitudes of the eigenmode of the 
largest eigenvectors on So, and (C) Fraction of +1 in Soo, as function of r. 



